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We study the effect of inhomogeneities in Hall conductivity on the nature of the Zero Resistance 
States seen in the microwave irradiated two-dimensional electron systems in weak perpendicular 
magnetic fields, and we show that time-dependent domain patterns may emerge in some situations. 
For an annular Corbino geometry, with an equilibrium charge density that varies linearly with 
radius, we find a time-periodic non-equilibrium solution, which might be detected by a charge 
sensor, such as an SET. For a model on a torus, in addition to static domain patterns seen at high 
and low values of the equilibrium charge inhomogeneity, we find that, in the intermediate regime, 
a variety of nonstationary states can also exist. We catalog the possibilities we have seen in our 
simulations. Within a particular phenomenological model, we show that linearizing the nonlinear 
charge continuity equation about a particularly simple domain wall configuration and analyzing the 
eigenmodes allows us to estimate the periods of the solutions to the full nonlinear equation. 



I. INTRODUCTION 

A novel "zero-resistance state" (ZRS) was observed 
a few years ago, when a two-dimensional electron gas 
(2DEG), subjected to a weak magnetic field, was also ir- 
radiated with microwaves! 1 * 2 * 3 * 4 * In the ZRS state, the dc 
resistance measured in a Hall bar experiment appeared to 
vaniish, while correspondingly, the longitudinal conduc- 
tivity a XXJ measured in a Corbino geometry appeared to 
be zero. A related phenomenon of microwave-induced re- 
sistance oscillations had been observed earlier in samples 
of lower quality, in which there are changes in the dc re- 
sistance produced by the microwave irradiation, which 
depend in an oscillatory manner on the ratio of the 
microwave frequency to the cyclotron frequency of the 
2DEG, but where the resistance was not driven down to 
zero.™ 

Following the discovery of the ZRS, a phenomenolog- 
ical explanation was put forward.^ It was assumed that 
for microwave power above a certain threshold, for an ap- 
propriate sample, in certain ranges of the magnetic field 
and microwave frequency, the microwave-induced resis- 
tance oscillations would increase to the point where the 
differential longitudinal conductivity would become neg- 
ative at zero dc electric field, and it was noted that in this 
case that the zero-electric-field solution would necessarily 
be unstable. It was argued that the sample would spon- 
taneously break up into a pattern of domains, separated 
by domain walls. Within each domain, the magnitude of 
the electric field would be a constant, E c , dependent on 
the strength of the microwave field and other parameters, 
which satisfies a condition that the longitudinal current 
generated by the field is zero. The direction of the electric 
field would change discontinuously at a domain wall, and 
might also vary continuously within a domain. It was as- 
sumed that the domain walls could move rather easily, so 
that, e.g., if a finite voltage is applied between the inner 
and outer edges of a Corbino sample, the system would 
respond by a displacement of domain walls, without caus- 
ing the magnitude of the field to deviate from the critical 
value E c , and without inducing a longitudinal current in 



the sample. 

If one accepts this general picture, an obvious ques- 
tion is what determines the domain pattern in any given 
sample? In fact, one may question whether there is nec- 
essarily a static domain pattern at all. Since the system 
is driven out of equilibrium by the constant absorption 
of microwave radiation, it is possible in principle that the 
system will enter a time-dependent state, where domain 
walls may move about in a chaotic, or possibly quasiperi- 
odic, manner. It is our purpose here to address this ques- 
tion. We shall argue that time-dependent states are likely 
to occur in at least some situations, and we shall propose 
one experimental geometry where it should be possible 
to observe this time-dependence. The dc properties of a 
time-dependent state should be similar to those of a state 
with static domain walls; in particular the observed value 
of a xx should be zero, or at least very small compared to 
the conductivity in the absence of microwave radiation. 

At least two different microscopic mechanisms have 
been proposed to explain the microwave-induced resis- 
tance oscillations at lower microwave intensities, and 
to produce the negative zero-field conductance assumed 
in the ma croscopic model of ZRS: a "displacement 
mechanism^ 8 ' 9 ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' in which the absorption 
of a microwave photon by an electron leads (for appropri- 
ate values of parameters) to a disorder-assisted displace- 
ment of the cyclotron orbit in an up-hill direction with 
regard to t he lo cal dc electric field; and a "population 
mechanism'^ 17 * 1 8 1, in which absorption of the microwave 
photons leads to a population inversion in the partially 
filled Landau levels close to the Fermi level. Although 
these mechanisms may lead to quite different estimates 
for such parameters as the threshold microwave power 
or for the critical field E c at a given level of microwave 
irradiation, it appears that they give rise to qualitatively 
similar phenomenological models. In any case, we shall 
not make any assumptions about the particular micro- 
scopic mechanism in this work. 

In formulating the equations of phenomenological 
model, it is convenient to divide the local electric cur- 
rent into a longitudinal or "dissipative" part, and a Hall 
current which is always perpendicular to the local field E. 
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In previous workp^D we considered a model where the 
functional form of the dissipative current could vary from 
one place to another in the sample, due to small inhomo- 
geneities in the doping density, or due to other sources 
of disorder. However, we assumed the Hall current to be 
governed by a linear (in field E) Hall conductivity, which 
we took to be uniform throughout the sample. In this 
case, we were able to derive a Lyapunov functional, whose 
value can only decrease in time, for a system with spec- 
ified electrochemical potential at the boundaries. Then 
the system must eventually reach a static state, whose 
potential configuration is at least a local minimum of the 
Lyapunov functional, so that oscillatory or chaotic time- 
dependent solutions are not possible as a steady state. 



If we allow the Hall conductivity to vary from one place 
to another, however, it is generally not possible to find a 
Lyapunov functional for the system. In this case, there 
is no minimization principle to determine the long-time 
behavior, and there is no guarantee that a stable time- 
independent solution exists. In fact, we shall explore ex- 
plicitly some simple geometries where, for certain ranges 
of parameters, the long-time behavior is time-dependent. 



The rest of the paper is organized as follows. In Sec. 
PH we discuss the phenomenological relation between the 
dissipative current and the local field, provide a brief 
overview of the Lyapunov functional formalism, then use 
the local capacitance model to arrive at a continuity 
equation that can be solved for the position-dependent 
electrochemical potential. We consider the solutions to 
the continuity equation for two geometries, the first be- 
ing Corbino-like, with the potential held fixed at two 
edges along one direction and periodic along the sec- 
ond direction. In a second geometry (torus), periodic 
conditions along both directions are used. Having estab- 
lished in Sec. [IT] that only variations in the equilibrium 
charge density contribute to the divergence of the Hall 
current, in Sec. [HI] we describe, for both geometries, the 
solutions to continuity equation for uniform equilibrium 
charge density. Sec. [IV] contains the central result of 
the paper; we show that, in a Corbino setup, linearly 
varying equilibrium charge density would give rise to pe- 
riodic solutions. In Sec. [VJ we present the numerical 
evidence for the presence of nonstationary solutions on 
the torus. The following section, Sec. [Vl| is dedicated to 
gaining analytical understanding of how solutions of pre- 
vious section arise by analyzing the linearized continuity 
equation. In Sec. |VII| we provide some details on how 
simulations have been carried out and compare the pre- 
dictions of the linearized analysis of the previous section 
with the numerical simulations. We comment on experi- 
mental prospects in our concluding section, Sec. |VIII| A 
discussion of some alternative assumptions for the Lya- 
punov function, and their consequences for the relative 
costs of different domain- wall configurations, is given in 
an Appendix. 



II. MODEL 

In the presence of an external microwave field, we as- 
sume the following relation between the dc current j(r) 
and the field E(r): 



j"(E,r). 



-a H z x E - AV 2 E, 



(1) 



where E = -VF(r) = V/ ^ r ) and z is a unit vector nor- 
mal to the plane. Here, p(r) is the electrochemical po- 
tential and — e is the electron charge. For the sake of con- 
venience, we will from here on refer to E as the electric 
field and V(r) as the electrochemical potential. In Eq. 
[I] we will explicitly allow the Hall conductivity <jh to be 
position-dependent. The dissipative current j d satisfies 
the condition that the differential dissipative conductiv- 
ity, 0-^0 (E) = djd/dEp, is symmetric: 



<3(E,r) = <(E,r) 



(2) 



The nonlocal (third) term in Eq. ([I]) implements an ul- 
traviolet cutoff, which will lead to a finite domain wall 
thickness, proportional to the square-root of the param- 
eter A, which will be taken to be small compared to the 
size of typical domain (indeed, this is almost a matter 
of definition of domains phase). In practice, we expect 
that the domain wall thickness will be of the order of 
the cyclotron radius l c , which is of order 1 pm in typical 
samples where ZRS is observed. 

The vector function } d may also depend explicitly on 
the position r, due to inhomogeneities in the 2DEG, and 
its direction may not be perfectly aligned with E. In 
previous work, we considered explicitly the effects of in- 
homogeneities in } d arising from gradients in the equi- 
librium electrostatic potential (j)^ due to disorder. In the 
present paper, however, we shall ignore this complication, 
and shall assume that the function } d has no explicit de- 
pendence on r, except for a brief discussion at the end. 

Eq. ([I]) will be supplemented by the continuity equa- 
tion, 



P : 



Vj, 



(3) 



where p is the charge density. Writing E = — VV(r), we 
may relate changes in the electrochemical potential V(r) 
to changes in p through the inverse capacitance matrix 
W: 



SV(r) 



J d?r'W(r,r')6p (r) . 



(4) 



If a time-independent steady state is reached, then we 
have simply V • j = 0, and the precise form of W is 
unimportant, but the form of W will be relevant for time- 
dependent solutions. 

In a Corbino geometry, one specifies the potential on 
the inner and outer boundaries of the sample, and one 
looks for a solution for V(r) consistent with these bound- 
ary conditions. If one assumes <jh to be a constant, the 
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Hall current cannot contribute to V • j in the interior of 
the sample, so it does not appear in Kirchoff 's equations. 
Consequently, the solution V(r) is independent of <jh and 
we may, for simplicity set cfh = 0. To recover the Hall 
current, one simply inserts the resulting solution for E 
into the second term in (JlJ, at the end of the calculation. 

Condition ([2J on } d allows us to define a scalar Lya- 
punov functional as 



G[V] 



d 2 r 

E(r) 



g(E) + ^ (V • E) 2 



g = / dE'-j d (E') 

JO 

A variation of ^ is given by 



SG= d 2 rV -j l 5V 



dsn -ySV, 



(5) 
(6) 

(7) 



bound 



where j l — } d — AV 2 E. The second integral vanishes on 
equipotential boundaries, or in the absence of external 
currents. Then, if cjh is independent of position, the ex- 
trema of G are found to be steady states, with V • j =0. 
Using the positivity of the inverse capacitance matrix W, 
one may show that G[V(t)] is indeed a Lyapunov func- 
tional, i.e. a non-increasing function of time, so that its 
minima are stable steady states. In general, G may have 
multiple minima. Any initial choice of V(r) will relax 
to some local minimum of G, though not necessarily the 
"ground state" with lowest G. Nevertheless, one might 
expect that in the presence of noise, the system might 
tend to escape from high-lying minima and wind up in a 
state with G close to the absolute minimum. 

We shall assume here that the function g depends only 
on the norm squared E 2 of the electric field, as would be 
appropriate in the limit of a uniform isotropic electron 
system, assuming that there is no unique axis picked out 
by external factors such as the polarization of the mi- 
crowave field. By hypothesis, under conditions where 
ZRS occurs, the function g must have its absolute min- 
imum at a non-zero value of the electric field, E = E c . 
Expanding about this minimum, and keeping only the 
nontrivial terms of lowest order, we may write 



5 (E) = g(E c ) 



8£2 



{E 2 - E 2 C 



(8) 



where the coefficient <j c has the dimensions of a conduc- 
tivity. In the absence of other information, we take <j c to 
be of the order of the dark conductivity (see Sec. |VIII 
for estimates). When E = E CJ the dissipative current, 
} d = dg/dE, will vanish, and cr^(E) = cr c E a Ep / E 2 . 

In the following calculations, for reasons of simplicity, 
and also to avoid introduction of additional parameters, 
we shall assume that Eq. ^ is exact for all values of E. 
(Actually, only the range < E < E c is important for 
our calculations.) We shall also assume that E c and a c 
are independent of position. 



The charge density p(r) is naturally broken into two 
parts, 



p(r) = p (r) + 5p(r), 



(9) 



where po(r) is the density in thermal equilibrium with 
no incident microwave radiation (by definition, the elec- 
trochemical potential is constant everywhere under such 
conditions). In a Corbino setup, for example, po(r) is the 
equilibrium density when the contacts at both edges are 
set to V = 0. In addition to having variations due to a 
non-uniform local dopant density, po( r ) could be tuned 
by a voltage on an external gate that is displaced from the 
2DEG by a distance which varies from one point to an- 
other. The remaining term, 5p(r), is the non-equilibrium 
charge density produced by microwave-induced domain 
structure, and it is this contribution which is responsi- 
ble for the position dependent electrochemical potential, 
according to Eq. Q . 

In our calculations, we shall assume a local form for the 
capacitance matrix W , so that the electrochemical poten- 
tial V(r) and the nonequilbrium charge density Sp(r) at 
a given point are simply proportional to each other: 



V(r) = C^Spir) 



(10) 



This form will be correct if there is a parallel conduct- 
ing gate, set back from the 2DEG by a distance d small 
compared to the typical domain size, in which case we 
have C « eeo/d, where e is the dielectric constant of the 
material between the 2DEG and the gate. If a nearby 
conducting plane is absent, the local capacitance model 
will not be strictly correct, but we would expect to ob- 
tain qualitatively correct results by using (10), replacing 



the setback distance d by a characteristic domain size. 

We shall assume that the Hall conductance at point r 
is determined in the usual way by the charge density p 
at that point : 



B 



Variations in the Hall conductance due to variations in 
po(r) will be crucial for the effects we investigate below. 
On the other hand, variations in cjh due to variations in 
Sp(r) play no role in the dynamics when one assumes a 
local capacitance, as in Eq. (10). The contribution of 



this term to is parallel to VV, and, therefore, gives 
no contribution to the divergence of the Hall current. 
Hence, we shall neglect this contribution and write the 
Hall conductance as 



Po( r ) 
B ' 



which is a quantity fixed at the outset, independent of 
the non-equilibrium charge-density induced by the mi- 
crowave radiation. 

We finally note that the Lyapunov functional of Eq. [5] 
misses one important piece of physics-while it incorpo- 
rates the fact that there is a high penalty for the field 
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to change o n leng th scales shorter than the domain wall 
thickness y/X/o^, it neglects the fact that, at least in 
the local capacitance model, the domains cannot be ar- 
bitrarily big. More precisely, if we imagine drawing lines 
through the domain that are parallel to the field inside of 
it, then these line segments inside the domain cannot get 
too long. Otherwise, the induced voltage difference be- 
tween points at opposite ends of these lines would result 
in an induced charge density dp that alters the total elec- 
tronic charge density by a large amount, which should 
be taken into account in the functional form of the Lya- 
punov function g(E). For example, one might need to 
take into account variations in the parameters E c and a c 
in Eq. (|8| due to variations in Sp. In a local capacitance 
model, the quantity 



-^max — q ' \- L - L J 

where po is the mean electronic charge density in the 
2DEG, will clearly provide an upper cutoff to the size 
of any domain. Rather than modifying the functional, 
we shall assume here that the domains we consider are 
below such cutoff. Note that the maximum domain size 
will tend to be largest if the setback distance d of the 
screening gate is large, so that the capacitance per unit 
area is small. 



Dimensionless parameters 

In our calculations below, we shall assume a rectangu- 
lar sample with linear dimensions L x = L and L y = kL. 
As postulated above, in the domains phase, the mag- 
nitude of electric field away from the domain walls be 
equal to E c . We will then work with dimensionless 
units, with distances measured in units of L, electric field 
E(r) = —W(y) measured in units of E c , electrochemical 
potential V(r) measured in units of E C L, conductivities 
measured in terms of the coefficient cr c , and, finally, time 
measured in units of CL 2 /a c . 

The dimensionless Lyapunov density, c.f. Eq. ([5|, we 
consider is then g(E) = (— ^- + ^-). When it is min- 
imized, the electric field is fixed at 1, that is E c , away 
from the domain walls. Including the nonlocal contri- 
bution giving nonzero thickness to domain walls, the di- 
mensionless longitudinal current density implied by the 
Lyapunov form we use is 

j'(25) = (£?-l)|-*f V 2 E, (12) 
where ld w is related to parameter A from Eq. [5] via = 

1 /a 

l y cr c * 

In the reduced units introduced above, we will be solv- 
ing ^ = — V • (} l + j^), for electrochemical potential 
V(r; t), on a [0, 1] x [0, k] rectangle. We have looked at 
two different sets of boundary conditions. In the first one, 



corresponding to Corbino geometry, we assume periodic 
boundary along x, but the potential is being fixed to be 
zero at the top (y = ft) and the bottom (y = 0) bound- 
aries. In the second case, corresponding to a torus, the 
periodic boundary conditions along both directions are 
assumed. For simplicity, will only consider unidirectional 
Hall conductivity <JH{y)- 



III. UNIFORM HALL CONDUCTANCE 

We first consider the case where the equilibrium charge 
density, and, therefore, Hall conductivity gh is uniform. 
Then the domain patterns observed should follow directly 
from Lyapunov energetics. If 7 is the angle that the elec- 
tric field forms with the domain wall, then, for the Lya- 
punov density that we have assumed, the cost of domain 
wall per unit length is proportional to sin 3 7 (See the 
Appendix for the proof). If we neglect the energy cost 
associated with wall crossings and exponentially small 
interactions between well-separated walls, then we can 
describe the ground states, that is the lowest Lyapunov 
energy states, of the system with the two sets of bound- 
ary conditions introduced above. They are illustrated in 
Fig. [1] 

In the simpler case of all-periodic boundary conditions, 
the ground state, shown in panel a) of Fig. [I] would have 
two horizontal and two vertical walls, with electric field 
forming angle 5 (tan 6 = k) with horizontal walls. Of 
course, there will be a degeneracy of these states due 
to translational invariance; the pattern can be displaced 
by a arbitrary amounts in the horizontal and vertical 
directions with no change in Lyapunov energy. 

The situation is more complicated in the case of the 
Corbino-like geometry, where we assume periodic bound- 
ary conditions in the x-direction, but the top and bottom 
boundaries are grounded. Results actually depend on the 
sample aspect ratio k. For k > 0.5, the case illustrated 
in panel b) of Fig. [I] the lowest energy state has a single 
horizontal domain wall at y = k/2. (There is actually a 
doublet of ground states, related to each other by a global 
flip of the signs of the electric field.) For k = 0.5, there 
are two degenerate ground states, the first again consist- 
ing of the single horizontal wall centered at the middle of 
the domain, with the second, shown in panel c) of Fig. 

t having four domain walls forming an angle of J with 
e top/bottom boundaries, and the electric field point- 
ing along ±y in the domains adjacent to the top and 
bottom boundaries, and along ±£ in the domains that 
are fully in the bulk. This second solution can, of course, 
be displaced in the horizontal direction by an arbitrary 
amount, without penalty. For the case of n < 0.5, there 
are an uncountable infinity of ground states, which may 
be thought of as the combination of the single horizontal 
wall and the diagonal wall states, as illustrated schemat- 
ically in panel d) of Fig. [I] The interactions between 
the walls will most likely break this continuum down to 
a finite set of states, but such investigation lies outside 
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Figure 1: Schematics of domain patterns in a rectangular 
sample with aspect ratio ^ = k for the case of uniform Hall 
conductivity. Thick lines are sample boundaries, with solid 
ones corresponding to boundaries where the electrochemi- 
cal potential is fixed, and the dashed ones corresponding to 
boundaries where periodic boundary conditinos are employed. 
Thin lines are domain walls, and the arrows indicate the direc- 
tion of the electric field within the domains. The magnitude 
of the field is constant away from the domain walls. In panel 
a), the periodic boundary conditions along both directions are 
assumed, and the shown state of the lowest Lyapunov energy 
has two horizontal and two vertical domain walls. Due to the 
periodic boundary conditions, this four-wall state is equal in 
energy to any state obtained by translating the original one 
along either x or y direction (or both). The rest of the pan- 
els correspond to the case of periodic boundary conditions 
along x, with both the top and the bottom boundaries being 
grounded. Panel b) corresponds to the case of k > .5, panel 
c) shows one of the degenerate states for k = .5, and panel d) 
represents one of a host of different possibilities arising in the 
case of k < .5. 

the scope of this paper. 

It should be emphasized that the existence of non- 
equivalent degenerate solutions in the case of n = 0.5, 
and the multiple degenerate solutions for n < 0.5 are 
a peculiar feature of the quartic Lyaponuv function (|8| 
that we have assumed in our analysis. In this case the 
Lyapunov energy per unit length of a domain wall with 
field-angle 7 = 7r/4 is smaller than that of a domain wall 
with 7 = 7r/2 by precisely a factor of 2y/2, just compen- 
sating for the difference in total domain- wall length in the 
two solutions. We believe that a more realistic assump- 
tion for the form of the Lyapunov function g would most 
likely favor the solution with diagonal 7r/4 walls. (See 
discussion in the Appendix.) Similarly, we expect that 
a more realistic Lyapunov function would favor a maxi- 
mum number of diagonal walls, and the shortest possible 
length of horizontal wall, for aspect ratios n < 0.5. For 
any reasonable choice of g, however, it is likely that mul- 
tiple metastable configurations will exist, for a range of 
values of ft, including configurations with a single hori- 



zontal wall and configurations with two or more diagonal 
walls. 



IV. CORBINO GEOMETRY WITH A 
GRADIENT IN a H . 

A particularly simple model with a nonuniform equilib- 
rium density (and, hence, nonuniform Hall conductivity) 
will illustrate the possibility of having a nonstationary 
state. Suppose that the Hall conductivity varies linearly 
with position, crniy) = oty. This, of course, makes sense 
in the context of a Corbino geometry, but not on a torus. 
With linearly varying cr^(y), we find 

divj^ = a' H (y) E x (r) = aE x (r). (13) 

Suppose that ^o( r ) satisfies divj d = 0. Then V(r,t) = 
Vo(r + atx) trivially satisfies ^ = — div(j d + J H ) = 
— olE x (y). In absolute units, this corresponds to a drift 
velocity 

v = ixVcr^/C. (14) 

Of course, this time dependence is not observable for a 
translationally invariant state, as in the case of a single 
horizontal wall in the center of the sample. However, if 
diagonal domain walls are present in the sample, then, if 
one were to detect electron density at a single point with, 
say, an SETpUH then one would find it changing with 
time. 

The question of whether the states other than the one 
with a single horizontal wall would be present in an actual 
sample is hard to answer definitively. Given a random ini- 
tial configuration for V(r), whether time evolution would 
take it to a state having diagonal walls is a matter of the 
relatives size of the domains of attraction for the such 
states compared with the domain of attraction for the 
state with one horizontal wall. Time-evolving the con- 
tinuity equation above numerically (more details on the 
numerics are provided later in the paper), we found that, 
for n < 0.5, starting with a random initial guess for the 
electrostatic potential, it is very likely that we end up in 
the state with diagonal walls. Moreover, as mentioned 
above, if we were to assume a more realistic form for 
the Lyapunov function g than the quartic form (|8|, we 
believe that solutions with diagonal domain walls would 
most likely have lower Lyapunov energies than the solu- 
tion with a single straight wall, in which case their basins 
of attractions would presumably be further enlarged. 

We may also consider the situation where there is a 
voltage difference V y between the top and bottom edges 
of the system. As long as V y is small compared to E c L yj 
the voltage can be accommodated by small displacements 
of the domain walls, giving a non-zero average of E y , 
without producing a net current in the ^/-direction. This 
can occur for the moving domains that we find when 
oh depends linearly on y in the same manner as for the 
stationary solutions appropriate to constant <jh- Thus, 
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the property of zero dc conductance, as measured in the 
Corbino geometry, is preserved for the time-dependent 
solutions. 



V. DYNAMICS ON A TORUS: NUMERICAL 
RESULTS 

While not as readily experimentally relevant as non- 
stationary states in a Corbino setup, it is nonetheless 
instructive to consider the kinds of nonstationary states 
we can get on a torus. We believe that the results might 
be applicable to cases where equilibrium density has cor- 
relation length that is much smaller than the sample size, 
and that, in such cases, one might get away with mod- 
eling a smaller piece of sample with periodic boundary 
conditions, rather than all of it with Corbino or Hall-bar 
boundary conditions. 

We consider here a sample of size L x — Ly — 1, with 
periodic boundary conditions. As described in Sec. |III| 
the equilibrium state for the case of uniform Hall conduc- 
tivity has two horizontal and two vertical domain walls 
with field-angle 7 = 7r/4. In our analysis, we shall assume 
that the spatially-varying part of the Hall conductivity 
has the simple form cr#(r) = asin27rn#?/, where n# is a 
positive integer. Our interests lie in trying to understand 
how the solutions depend on the non-uniformity param- 
eter a, for different uh and domain wall thicknesses 1& W . 

As we tune nonuniformity a away from 0, the Hall 
current develops a nonzero divergence in the equilibrium 
state. We find that, so long as the a is small enough, 
the four-wall solutions can accommodate such nonuni- 
formity by bending. We may imagine that the gradient 
of the <jh leads to a horizontal "drag force", proportional 
to z x V<7#, which can be counteracted, for small a, by 
the "restoring force" produced by a distortion of the do- 
main structure, away from the shape that minimizes the 
Lyapunov functional. By contrast, for large values of a, 
we find only solutions which are translationally invariant 
along the x-direction, and have have two or more hor- 
izontal domain walls with 7 = tt/2. The divergence of 
Hall current trivially vanishes for such states. 

What happens for intermediate values of a depends on 
Tin and 1& W . We find that if 1& W is big enough, then one 
can transition from the four- wall solution to the two- wall 
one without ever encountering the nonstationary solu- 
tions. As seen in Fig|2j the potential contour lines that 
have four-fold symmetry at a = stretch and rotate as 
a is tuned up, until they become horizontal. 

We find that increasing nn with everything else held 
fixed similarly suppresses the appearance of nonstation- 
ary states. At ldw ^ -H, however, nonstationary states 
do exist. For example, when l^ w = .028 and n# = 1, we 
find nonstationary states for a G [.07, .7]. Oftentimes, 
stationary and nonstationary states (perhaps even of dif- 
ferent types, see below) can exist at the same value of a. 
Which one we find in our simulations is both a matter of 
the initial conditions, as well as the relative size of the 
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Figure 2: Numerically-obtained electochemical potential 
V(r) for domain wall thickness ldw = -11, tih = 1. All V(r) 
shown are stationary. The transition between four-wall so- 
lution at a = and the one with two walls for a — 0.3 is 
continuous; that is, it is not interrupted by the appearance of 
nonstationary solutions at the intermediate values of a. 



basin of attraction for that particular solution. 

The nonstationary states we have seen can be broken 
into three groups. The first is intimately connected to 
the solution one gets at small a. The periodic solutions 
of the first kind can be thought as the system tunnel- 
ing from one small a state to the one connected to it by 
translation by \x. Figure shows snapshots of such a so- 
lution, plotting V(r, ti) for ti = i^, with i = {0, 1, 7} 
and T being the period, so that V(r, t + T) = V(r, t). 

In the periodic solutions of the second kind, the time- 
dependent potential is confined to only half of the sam- 
ple, and has the form V(r — xvt), representing a structure 
which drifts to the right or to the left with a velocity v. 
This behavior is, of course, similar to what we found in 
the Corbino geometry with a uniform gradient of an- 
Indeed, the observed values of v are, to a good approx- 
imation, related to the average value of Va# across a 
domain in the same way as the drift velocity was related 



to the gradient of cfh in Eq. (14) of Sec. IV 



Finally, the periodic solutions of the third type, seen 
only for n# = 1, can be thought as the combination 
of two solutions of type II: in one half of the sample, 
the potential moves to the right, while in the other half 
it moves to the left, with more complicated processes 
happening at the boundary between the two regions. 

As examples of these behaviors (all but the first one 
not shown in the figures), for a = .1, ldw = -028, we 
have found periodic solutions for uh = 1, 2, 3, with the 
periods being 2.49 (n# = 1, type I), 1.90 (n# = 2, type 
II), 1.64 (n H = 3, type I), and 1.21 (n H = 3, type II). 




Figure 3: Numerically obtained electrochemical potential V(r, t) for domain wall thickness ldw = .028, tih = 1, a = 0.1 at 
different points during its period. The solution here "tunnels" from one a = solution to the one related to it by translation 
by x/2. 



VI. ANALYTIC RESULTS ON A TORUS 



A. Linearized Continuity Equation 



We would like to understand what controls the periods 
mentioned towards the end of the previous section. The 
quartic Lyapunov density of Eq. [8] is assumed through- 
out this section. It turns out to be useful to analyze the 
stability of the simplest solution, the one with two do- 
main walls perpendicular to y-axis. With the walls fixed 



at y 



(we keep the aspect ratio k at 1), and the 



field varying along y-axis only, the continuity equation 
is satisfied for all values of nonuniformity a, as well as 
all values of n# and domain- wall thickness ldw However, 
the numerical evidence cited above at the very least seems 
to suggest that, as we decrease a, the basin of attraction 
for the two- wall solution rapidly shrinks. We will demon- 
strate below that, below some a c (ld w ), the two wall solu- 
tion becomes unstable. As we have no analytic solution 
to the full-blown nonlinear partial-differential equation, 
the technique of choice will be to linearize the differential 
equation satisfied by the potential around the two-wall 
equilibrium and see whether the eigenmodes grow or de- 
cay with time for various values of a and domain wall 
thickness ldw After linearizing the equation, we will sep- 
arately consider the effects that nonzero ld w and nonzero 
a have on the nature of eigenmodes and the correspond- 
ing eigenvalues. We will then try to understand what 
happens when both of the parameters are nonzero. 



Linearizing around the two- wall solution, we see that 
changing the field by 5~E changes the longitudinal current 
by 



El 



1 



SE + El- 5E V 



J 2 



Here, Eo(y) is the field corresponding to the solution with 
two horizontal walls. For the quartic Lyapunov density, 
the form of E (y) is derved in the Appendix for the case 
of a single domain wall. That result should be applica- 
ble in the current case of two domain walls, so long as 
the separation between them greatly exceed domain wall 
thickness, ld w <C .5. The change in the Hall current is 
6} H = a H (SE y x + SE x y). Employing <SE = -V5V, we 
see that the corresponding divergences turn out to be: 



divJj 



-dydVdyEl 



EidyySV+JfV^SV. 



Furthermore, since gh depends only on y, we find 
divS} H = (d x 5V)(d y (jH)- The charge continuity equa- 
tion, aided by the assumption of the local capacitance 
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model, then tells us that 

d t 5V = -divSj = ] ^t-±V 2 5V + I (dySV) d y E 2 

+e 2 dy y sv - l ^fvHv + (a y a H ) a^v. 

Different Fourier modes along x-axis will be decoupled, 
motivating the ansatz 5V(x, y) = e at e 27vnxt h(y). Our dif- 
ferential equation then becomes: 



Aj2 

dw 



il^Y^-M^h-ld.Elh' 



3^ 2 " 1 



o 2 / 2 , 

1x''dw \ foil 



(15) 



where q x = Inn. 



boundary layers are most conveniently written as func- 
tion of z = (y—y W aii)/ldwj the matching conditions would 
be 



Ubulkf + 

dh\ ulk (y) 



wall / 



/if dn/ (oo) 



iibndry 
an i+l 



y y-waii 



dz 



etc. 



Assuming the quartic form (Eq. [8| for the Lyapunov 
function g(E), we have worked out the first few lead- 
ing terms in the expansion of a sym , the eigenvalue cor- 
responding to the mode that is invariant under reflec- 
tion through both domains wall and domain midpoints 
(that is, h$ ornainl = h$ omain2 ) and a asym , the eigenvalue 
corresponding to the mode that is invariant under re- 
flection through domain midpoints but acquires a minus 
sign when reflected though the domain walls (that is, 
h domaini = _ h domain2y Ag the com p uta tion is rela- 
tively straightforward, we will omit the details of it here 
and present the final result: 



B. Case of a = 

We first consider the case of uniform Hall conductivity 
but nonzero domain- wall thickness. Equation (15) now 
turns into 



aH 2 

Qx\ 



x dw 



3E 2 



Qx 



2 E 



a 2 l 2 



x l dw 



h" 



-dyE 2 h' 

l lw h nn = Q 



The regions immediately surrounding the domain wall 
and the regions far away from it should be considered 
separately, turning this problem into a WKB-type com- 
putation. The interval [0, 1] is broken into four regions, 
two bulk domain regions and two thin (extending a few 
Idw) regions around the domain walls at \ and \. The 
equation is separately solved in each region, then differ- 
ent solutions are patched at various boundaries. We will 
assume a power series expansion for eigenvalue a and 
eigenfunction h(y): a = Y. a i l d w i %) = T, h i{y) l dw 
Trivially, we have ao = and ho = C. 

Note that E 2 is symmetric under reflection through 
either the center of domain wall or the center of the 
domain. That means that we can pick the solution to 
be symmetric or antisymmetric with respect to reflection 
through the centers of domains or through the domain 
walls. There are four cases to consider, and we will first 
briefly mention those two where the solution is antisym- 
metric under the reflection through domains' midpoints. 
Going order by order in the expansion, one can show that 
all contributions are identically zero. We now turn to the 
more interesting case of midpoint-symmetric solutions. 

Solving for h^s in the bulk and in the boundary lay- 
ers, and matching these solutions yields the terms in 
the eigenvalues expansion. We note that since hi in the 



V/ 2 

^Vx^dw 



dw 



(16) 



These expressions illustrate two points. As expected, 
for high-g^ modes, both eigenvalues are negative, simply 
reflecting the fact that modulation of potential on wave- 
lengths smaller than domain wall thickness are strongly 
suppressed. For relatively small q x , though, a sym > 0. 
In particular, considering q x = 27r, time evolution of the 
eigenmode (see lower right panel of Fig. [5]) takes us from 
the potential that has two horizontal walls to the solution 
that we know to be stable, the one with two horizontal 
and two vertical walls, i.e. panel a) of Fig. [l] 

Finally, we emphasize that the periodic conditions 
along y caused the zeroth order (in q 2 ldw) term to van- 
ish. In a Corbino setup, the zeroth order term is neg- 
ative, thus causing both eigenvalues to be negative for 
thin enough domain walls. This implies that the solution 
with two horizontal walls would be locally stable in the 
Corbino geometry. 



Case of 1 



In the limit of zero domain wall thickness, E 2 



1 ev- 
erywhere but on the domain walls, where O v Eq diverges. 
Setting Idw = in Eq. [15] must then be accompanied by 
the prescription of how to deal with h(y) on the domain 
walls. This is not a straightforward task. Instead, we 
shall take a stab by assuming a simple boundary condi- 
tion, namely that h! — as one approaches the domain 
walls. This boundary condition will suffice to determine 
the eigenvalues, as (15) is now simplified, away from the 



domain walls, to a second order differential equation, 
( a — 2i\q x OL% cos 2n (y + (3) ) • h = h" . 
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(Since the positions of domain walls are fixed, we have 
introduced an offset parameter f3 to study the general 
case where the domain wall is not necessarily located at 
an extremum of cfh-) 

We shall be interested in understanding the behavior 
of eigenvalue a(a) for small a; in particular, we will be in- 
terested in the most physically relevant eigenvalue, with 
the largest (most positive) real part, as a — > 0, a(a) — > 0. 
Note that the requirement that h! vanishes at each do- 
main wall essentially decouples different domains. The 
goal is to work out the expansion of eigenvalue a and 
eigenfunction h(y) in powers of a: a = ^aja^ and 
h(y) = ^jhj(y)ai . One readily sees that real part of 
eigenvalue (and eigenfunction) only contains even pow- 
ers of a, while the imaginary part has only odd powers. 

Starting with clq = and ho(y) = 1, we can get h\(y) 
by integrating twice the differential equation it obeys, 

a\ — 2irq x cos 2ir{y + /3) = h". 

The assumed boundary condition, h[ = 0, then 

fixes a\. In a similar manner, and higher order 

terms in the expansion can be obtained. The first few 
terms in the eigenvalues expansion are as follows: 



D. 



Case of Idw, ol ^ 



ai 

a 3 



i4q x cos2tt(3; 
4(2-% 



12 * 

2 ) cos 



ql cos 



45 



2irp ( 1 



120 

^2" 



1080 



/1080 435 1 



1 



cos 4tt(3 . 



The next term in the expansion would not be given ex- 
plicitly here, but we find it useful to give numerical values 
of the a[s here: in the case of /3 = 0, n = 1, one finds 
a t = 25.13, a 2 = -.8987, a 3 = .0149, a 4 = -.0041. We 
also find that the first two terms in the expansion of real 
part of a, d2(P)a 2 + a^fya 4 < for all (3 and for a 
smaller than about 8.36. 

We finally comment that making shift (3 — > (3 + \ is 
equivalent to moving from one domain to the other. This 
leaves the real part of eigenvalue unchanged while flipping 
the sign of the imaginary part. It follows that eigenvalues 
come in conjugate pairs, and that corresponding eigen- 
functions must identically vanish in one of the domains. 
This last result, of course, is a consequence of our as- 
sumption that h! = at the domain walls, but it is remi- 
niscent of the type III periodic solutions in our numerical 
results, where charge is moving in one directions in one 
of the domains, and in the opposite direction in the other 
domain. In this light, it is not at all surprising that the 
solutions of this kind show up at fairly large a, where the 
fact that the domain wall thickness is actually nonzero 
is not qualitatively important. This gives us reassurance 
that setting h'(y) = on domain walls was a reasonable 
guess. 



We are now ready to address the case where the Hall 
conductivity is nonuniform and the domain walls have 
nonzero thickness. We restrict the analysis to the cases 
where a is at most 0(1) and we consider the mode with 
the lowest q Xj q x = 2tt. The analysis of the previous 
subsection then implies that, for a/ld w ^> 1, we might 
expect 5ft a(a, ld w ) < 0. The eigenmode then is a decaying 
one, and we find that the solution with two horizontal 
walls is stable. Below some a c (which is a function of nn 
and ldw)i however, the mode will start growing, and we 
find that putting (just about any) perturbation on top 
of the two-wall solution and time-evolving the resulting 
guess would take us away from it. Indeed, for a/ldw <^ 1, 



\Uw) > (see Eq. 16) 



we expect a(a,ld w ) 

There is, additionally, another special value of a, a^, 
below which the eigenvalues are real and unequal to each 
other, whereas above it the eigenvalues are complex and 
are conjugate to each other. At that special point, the 
eigenvalues are, of course, equal to each other. We ex- 
pect that Oidjldw is 0(1). We, however, can't estimate 
fl(^dj'diu); thus we can't establish whether a c > ad or 
the other way around. It is plausible, though, that non- 
stationary solutions would exist for a G [a^,a c ], as, in 
that case, we would have oscillatory runaway modes. 

As for the functional form of the eigenmodes, we expect 
them to be mostly symmetric between two domains for 
a < ad- If our guess for the boundary conditions in 
the previous subsection were a reasonable one, we would 
expect the modes to have support in one of the domains 
only for a/ldw ^> 1- 



VII. SIMULATION DETAILS 

A. Numerical Setup & Methodology 

The continuity equation was discretized on a triangu- 
lar lattice of grid points and then evolved in time. For 
the linear analysis part, the initial guess consisted of a 
two-wall solution with perturbation that only had one 
Fourier component in the direction parallel to the walls 
(along x). Though the precise form of y-dependence of 
the perturbation did not matter much, we have tried, 
among other things, a function that is constant on both 
domains and one that vanishes in one domain but is con- 
stant in the other one. In Fig. |4j we have plotted the 

' 2 



dpi 
dt 



with 



root-mean-square amplitude A(t) = 

the sum being taken over all grid points, as a function 
of t, for a particular solution. As can be seen in the 
plot, for small values of t the amplitude increases ex- 
ponentially over time, while also exhibiting oscillations. 
From this portion of the plot, we can infer both the real 
and imaginary parts of the eigenvalue a/, which governs 
behavior in the linear regime. At larger times, the ex- 
ponential growth has saturated, but oscillations persist. 
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Figure 4: Time trace of A = y ^ . | ^ | , with the sum being 

taken over grid points, can be used to extract both the real 
and imaginary parts of frequencies. The transition from linear 
regime (where the real part of the frequency co\ is nonzero) to 
the saturated one (where the real part of cj sa t vanishes) is 
apparent. 



a =1.0 



a =.7 




a =.3 



a =0 



in 



From this portion of the plot, we can identify a second 
frequency, a sat which characterizes oscillations in the sat- 
urated regime, as indicated in the figure. 

When we were looking for various possible nonlinear 
solutions, the initial guesses corresponded to the poten- 
tial at each site being assigned a random number between 
-0.5 and 0.5. Most of the analysis was carried out on a 
lattice with mesh size though resolution studies have 
also been carried out on grids with mesh about half as big 
and quarter as big. Working with the rougher mesh has 
restricted us to domain wall thicknesses ldw no smaller 
than about .03. 



B. Linear Regime Results 

We first comment on the solutions and then on the 
eigenvalues. Recall that the prediction was for eigen- 
modes to be largely localized to one domain for large 
a and be symmetric between the two domains around 
a = 0. And this is indeed what is seen in Fig. [5] The 
presence of free constants in the perturbative expansions 
we obtained for the eigenmodes prevents us from making 
more detailed comparison between theory and simula- 
tion. 

The eigenvalues inferred from numerical simulations 
in the case of ldw = -028, nn = 1 are shown in Fig. [6] 
For a = case, we get an excellent quantitative agree- 
ment between the simulation results and the contribu- 
tions of the first few terms in perturbative expansion for 
the eigenvalue a, shown in Eq (16). As predicted, the 



eigenvalue remains real for a while before turning com- 
plex, and these complex ones do indeed come in pairs. 
For this particular domain wall thickness, we find, using 
the notation of Sec. |VID[ ad ~ 0.07, a c w 0.7. We have 

is sufficiently big, the order is 
Note that the cur- 



Figure 5: Plots of absolute value of eigenmodes of the lin- 
earized continuity equation for all-periodic boundary condi- 
tions, ldw = .028, and nn — 1. The eigenmodes are devia- 
tions from the solution with two horizontal walls; thin dashed 
lines indicate the location of these domain walls. Plus and 
minus signs indicate whether the deviation is positive or neg- 
ative, while white color corresponds to zero deviation. For 
sufficiently big values of the nonuniformity a, all of the eigen- 
mode's support is within one domain. On the other hand, 
for a = 0, the (unstable) mode is symmetric between two 
domains. 



found though, that if 1& W 
reversed and one would have otd > a c 



vature of a(a) for a > ad is quite close to a2 of Eq. 17 
(solid line of Fig. [5). We can't predict the vertical offset 
between the simulation results and the solid line-this is 
a direct consequence of our inability to predict ad and 
a(ad) which we have referred to earlier. 

As for the imaginary part of the eigenvalue, the agree- 
ment between the simulation results and the perturbative 
expansion is excellent. It is noteworthy that not only is 
the agreement excellent in the linear regime but even by 
the time we get to the saturated regime, the agreement 
is still very good. Indeed, it was the this striking linear 
(in a) behavior exhibited by the full nonlinear solutions 
that has motivated the whole expansion enterprise. 



VIII. CONCLUSIONS AND EXPERIMENTAL 
PROSPECTS 

We have seen that non-uniformities in the Hall 
conductance <t#, produced by inhomogeneities in the 
equilibrium electron density, can lead, under certain 
circumstances, to time-dependent domain patterns in 
the microwave-induced Zero-Resistance states seen in 
2DEGs. Inhomogeneities in the equilibrium density may 
result from disorder, such variations in the density of 
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2.02, the sum of 
first two terms in 
the expansion 
in oc=0 case 
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Figure 6: Real and imaginary parts of eigenvalue a. The left panel shows the real part of a extracted from linear regime 
simulations on rougher mesh (filled circles), finer mesh (crosses), as well as the first term in perturbative expansion for the 
case of Idw = (solid line). The dashed lines connecting crosses and filled circles are merely a guide for the eye. Right panel 
shows the imaginary part extracted from simulations in the linear regime (with filled circles corresponding to rougher mesh 
and crosses to finer mesh), the nonlinear regime (open circles), and the first term in the perturbative expansion for the case of 
Idw = (solid line). 



ionized Si donors in the set-back doping layer, or from 
density variations that are deliberately imposed on the 
sample. Density variations of the latter type could be cre- 
ated by applying voltages to one or more gates, located 
at appropriate distances from the sample. 

In our calculations, we have studied explicitly two sim- 
ple models: a Corbino-type geometry, with a linear de- 
pendence of (Jh on the distance y form one edge of the 
sample; and a model with periodic boundary conditions 
and a variation in cfh that has a sinusoidal dependence 
on one of the coordinates. In both cases, we have found 
time-dependent solutions for appropriate choices of the 
parameters. 

The model with periodic boundary conditions may be 
useful as a first step to examine the effects of random 
density fluctuations in a macroscopic sample. In this 
case, we might identify the size L of our model with a 
region who size is given by the correlation length of the 
most relevant density fluctuations in the physical sam- 
ple. In our model we found that time-periodic solutions 
occurred for density fluctuations falling within a certain 
range. We speculate that similar time- varying solutions 
may occur in the random system, with some regions of 
the sample trying to generate time-dependent solutions of 
varying frequencies, and other regions favoring decay to a 
time-independent solution. If a time-dependent solution 
occurs in a macroscopic disordered sample, presumably 
the power spectrum of the dynamic potential fluctuations 
will contain many frequencies, and the domain- wall mo- 
tion will be chaotic. Since our starting equations are only 
valid on a length scale greater than the cyclotron radius, 
of order 1 /im, we would not want to consider intrinsic 
density fluctuations on the scale of the set-back distance, 



but would restrict ourselves to longer wavelength fluctu- 
ations of an extrinsic origin. 

The Corbino model studied in Section [IV] suggests a 
way in which time-dependence may by implemented in 
a controlled way, which may be most promising for ex- 
perimental realization. Experimentally, one might aim 
to create a circularly symmetric Corbino annulus, with 
an electron density gradient along the radial direction. 
The density gradient could be be produced by applying 
a voltage Vq to a front gate, whose distance s from the 
2DEG varied linearly with the radial coordinate r. If e s is 
the dielectric constant of the spacer between the gate and 
the 2DEG, then the imposed variation in the equilibrium 
electron density would be given by 



Po(r) - po 



s(r) 



(17) 



where po is the mean charge density at Vq = 0. For ex- 
ample, if s were to vary between 1 /im and 2 /im, from 
the inner radius to the outer radius of the annulus, and 
if e s w 13, one would need a gate voltage Vq = 0.09 V 
to produce a density difference of 3.5 x 10 9 cm -2 across 
the sample, which would be 1% of the average density 
for the sample studied in Ref. [2j The gate, as well as the 
spacer layer, should be transparent to the incident mi- 
crowaves; this might be achieved, for example, if the gate 
is made from a thin film with high resistivity. Alterna- 
tively, the density gradient in a 2DEG could be created, 
even if the spacer thickness is uniform by depositing sev- 
eral separately-contacted gates on top of the spacer and 
putting different voltages on them. 

We can use the above numbers to estimate the charac- 
teristic frequency of the charge oscillations that might be 
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observed. Let us consider a Corbino sample, whose width 
w is much smaller than its circumference L, say .1 mm 
and 10 mm respectively. In this regime, the curvature 
of the annulus is not important, and the periodic states 
should resemble those described in Sec. [IVJ To apply the 
calculations of Sec. |IV} we estimate the capacitance per 
unit area C = ee^/dhy equating the distance d to 1.5 
jam, the average distance to the gate, and use e s = 13 
for the dielectric constant. If we assume that the elec- 
tron density changes by a factor 1 + (3 across the width 
w of the annulus, then the drift velocity v in Eq. (14) 
has magnitude v = /3(j xy /wC. If we assume that the do- 
mains have a typical length of order w, then the moving 
domain structure will lead to a time-dependence of the 
charge density, at any given point, with a characteristic 
period given by 



Cw 2 



w 



(18) 



For the parameters considered here, if (3 corresponds to 
a density change of 1%, we find r p w 15 ns. 

As we have assumed that the domain size is set by 
the width w of the Corbino sample, to be consistent we 
should check that w = 100 /im does not exceed the cut- 
off scale L max = p/CE c of Eq. 
estimates of E c 

typical experimental density n = 3 x 10 n cra~ 
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The re ar e a few 
in the theoretical literatur e! 8 1 18 !; for the 

" 2 , the high- 
est estimate would be on the order 20 V/cm. Experimen- 
tally, while E c has not been measured, some of the results 
of Refs. I3|4I can be interpreted to imply E c ~0.1 V/cm. 
Cautiously, using the higher estimate for the critical field 
and the same estimate for capacitance per unit area C 
as above, we arrive at L max ^0.6 cm. So, in a Corbino 
sample that is 100 /im wide, the size of the domain should 
indeed be set by the the Corbino ring's width. By vary- 
ing the sample width, the spacer thickness, or the density 
gradient the time scale of interest, r pj can be tuned over 
a wide range. 

For the sample size we are considering here, we need 
to detect density changes at frequency of about 70 MHz. 
This might be accomplished in number of ways. An ap- 
porpriately tuned RF-SET, positioned on the surface of 
the the sample, underneath the gate should be able to 
pick up the fluctuations at that frequency. Alternatively, 
a small capacitor, placed at the surface, might be able 
to pick up changes in density. Finally, if one could put 
down the Ohmic contacts in the middlle of the sample 
without destroying Zero-Resistance state, then monitor- 
ing the voltage difference between the inside of the sample 
and its boundary, held at constant potential, might show 
the time dependence we are predicting. The magnitude 
of voltage fluctuations produced by the domain of width 
w is given by E c w. For a domain of width w — .1mm 
and E c = .lV/cm, this would imply voltage varying by 
lmV. 

The parameter a c in the Lyapunov function g(E) does 
not enter directly in our estimate for the oscillation 
period in the Corbino geometry. However, it should 



play a role in a system with disorder, as it did in our 
model with periodic boundary conditions, where the non- 
uniform drift velocity v is in competition with the ten- 
dency of the Lyapunov function to favor an incompat- 
ible domain structure. The dimensionless parameter a, 
which controlled the transition from time-independent to 
time-dependent solutions in Sec. [VJ may be written as 
a = /3cr X y/<j c , where (3 is here the fractional change in 
the electron density between its minimum and maximum 
values. The value of a c depends on details of the micro- 
scopic model, but it is expected generally to be of the 
order of the longitudinal conductivity cr xx in the absence 
of microwave radiation. For the sample described in Ref. 
12 in a field of ~.11T, where the most prominent ZRS 
was seen once the sample was subjected to microwaves, 
we estimate a value for cr xx of order 5 x 10 _5 S, while 



ff xy «5x 10~ 3 S. 



Finally, we note that our model for the Corbino geom- 
etry includes an intentional density gradient, but ignores 
effects of random density variations that may occur in 
an actual sample. If the random variations are too large, 
it is possible that they would pin the domain walls at 
specific locations in the sample, and defeat the tendency 
for motion produced by the intentional gradient. At the 
present time, we have little knowledge about the magni- 
tude of density variations over the range of length scales 
of interest and we cannot reliably estimate the impor- 
tance of any such pinning effects. It seems likely that 
pinning effects could be overcome by employing inten- 
tional density gradients that are sufficiently large. 
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Appendix A: ALTERNATIVE FORMS FOR THE 
LYAPUNOV FUNCTION. 



The Lyapunov energy cost of a domain wall depends 
on the form of the Lyapunov function #(E), for E < E c , 
as well as on the angle 7 between the domain wall and 
the field E just inside the domain on each side. (The 
angle must be the same on both sides of the wall.) For 
the maximum angle 7 = 7r/2, the field inside the domain 
wall passes through zero at the center of the domain wall, 
so that the entire range < E < E c is reached within the 
wall. If the field is not perpendicular to the wall, however, 
the component parallel to the wall remains constant and 
non-zero throughout the wall, and only the perpendicular 
component of E changes sign. The domain wall is then 
determined by the form of g in the range E c cos 7 < E < 
E c . 

In our calculations, we have used Lyapunov energy 
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density (|8|, which is a quartic function of the field 
strength E. It is the simplest possible function that has 
the following desired properties: it is an analytic func- 
tion of E 2 , and hence an an isotropic analytic function 
of E; it has an absolute minimum at a non-zero value 
of E, and a local maximum at E — 0. The function is 
then determined by the position E c of the minimum and 
the curvature a c at E = E c . However, if the quartic 
function is extended to values of E much larger than E c , 
the longitudinal differential conductivity, d 2 g/dE 2 is seen 
to increase without limit, going eventually as E 2 , which 
does not seem reasonable. It might be more reasonable to 
assume that at large values of E, the longitudinal differ- 
ential conductivity should approach a constant value a^. 
Although we are not directly interested in the behavior of 
g(E) for E > E c , the assumption that g is analytic does 
introduce a link between this behavior and the behavior 
at E < E c . In particular, if we assume a finite value 
for (Too, this implies that d 2 g/d(E 2 ) 2 — > for large E 2 , 
while it is positive at E = E c , which at least suggests 
that d 3 g/d(E 2 ) 3 is likely to be negative at E = E c . If we 
assume that g has the form of a cubic polynomial in E 2 
for E < E c , this would imply that g(0) — g(E c ) is larger 
than in case of a pure quartic form, for the same values 
of E c and a c . 

In previous work, (Refs. 19 20), we assumed the fol- 
lowing form for g: 



the longitudinal current j z is given in Eq. (12). 

Since |E| = E c at y = ±oo, the current density van- 
ishes far away from the domain wall. It follows that j l 
must vanish for all y. Then, the equation to solve is 

(Tf ^ ~ Ey AV Ey 



2E 2 ~ y 
?-£ c 2 sin 2 7 



2E 2 



E y - AV 2 'E y = 0. (A2) 



This is one of the few nonlinear differential equations that 
happens to be exactly solvable. The solution that satis- 
fies both the boundary conditions and the requirement 
that the domain wall be centered around y = is 



/ sin t \ 

E y (y) = £ c sin 7 tanh \-^=y\ (A3) 

Using this field profile, we find the Lyapunov cost per 
unit length of domain wall is given by 



JG( 7 ) 



[g(E) + ^-E) 2 -g(E c )]dy 



2^/\o~ c E 2 c sin 3 7 



(A4) 



9P (E) 



1 



(Tr.E: 



(1 + /3 2 



In 1 



(Al) 

The high field conductivity is then given by = <r c (l + 
f3 2 )/2. In the limit (3 — > oo, the function gp reduces to 
the quartic form g quarticOfy given by We note that 
(3 = 0.1 was used in numerical calculations whose results 
are shown in Figs. 1 and 2 of Ref. [19] and Figs. 4 and 5 
of Ref. M 

We shall now examine the way in which several dif- 
ferent assumptions for g affect the relative Lypunov cost 
of domain walls of different angles 7. We first demon- 
strate that, for the quartic form used in our calculations, 
the Lyapunov energy of a domain wall is proportional 
to sin 3 7. We will consider an infinite system, with a 
single domain wall centered around y = 0. The prob- 
lem is invariant under translations along x-axis, so that 
electric field E, and, consequently, longitudinal current 
j z , are functions of y-coordinate only. Moreover, the 
x-component of field is constant: E x = E c cos^. The 
boundary conditions on the y-component are £^(±00) = 
±£? c sin7. We are looking for a time-independent solu- 
tion of a continuity equation. Since both electric field 
E and Hall conductivity cjh depend on y only, the Hall 
current has a vanishing divergence. Therefore, solving 
divj = would only entail solving dj l y (y)/dy = 0, where 



A closed form expression cannot be obtained for the 
Lyapunov density given by Eq. |A1| In that case, we 
can solve j l y (y)= numerically, then integrate to get 
SG(j). Figure [7] shows the Lyapunov cost per unit length 
of a domain wall as a function of angle 7, normalized 
by its value at 7 = 7r/2, for quartic Lyapunov density 
9quartic(E>), as well as for 0/3- i(E) (used in Refs. 119120)) 
and gp=i (E), which gives = a c . 

Two points need to be made. First, we note that the 
ratios 6G(Tr/4)/6G(ir/2), for both (3 = 0.1 and (3 = 1, 
are smaller than the value 2 -3 / 2 that is obtained from 
the quartic form. Thus as stated in Sec. |III[ Fig. [7] 
demonstrates that more realistic forms of Lyapunov den- 
sity would break the degeneracy between having a single 
horizontal domain wall and a couple of diagonal ones, 
with diagonal walls (corresponding to 7r/4 domains) be- 
ing preferred at least in these examples. That, of course, 
is helpful, because it is the time evolution of 7r/4 walls 
that produces dynamics. 

The second point is that the energy cost of very small 
angle domain walls is actually universal, independent of 
the the precise form of g(E), if E c and a c are specified. 
Indeed, if the angle 7 <C 1, then even the smallest value 
that (E/E c ) 2 attains, cos 2 7, is very close to 1. But, 
by construction, all g(E) look alike in the neighborhood 
of E = E c . Hence, the ener gy penalty for a small angle 
domain wall is given by ( A4), irrespective of form of g(E). 
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Figure 7: Normalized Lyapunov cost of domain wall per unit 
length, SG(j)/SG(7r/2), as a function of the angle 7 between 
the field in the bulk of the domain and the domain wall. Solid 
line is the exact result for quartic Lyapunov density. Triangles 
correspond to the Lyapunov density gp of Eq. |A1| with the 
value (3 = 0.1, used in our previous work (Refs. 119120)) . while 
crosses show results for /3 = 1. 
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